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Abstract- In this paper we develop an analytic expression for calculation of the gradient 
of the acoustic pressure from a rotating blade on a moving surface for application to the 
Fast Scattering Code (FSC). The analytic result is intended to be used in the helicopter 
noise prediction code PSU-WOPWOP. One of the goals of the derivation is obtaining a 
result that will not use any more information than are needed for the prediction of the 
thickness and loading noise. The result derived here achieves this goal and its incorpora- 
tion in PSU-WOPWOP is straight forward and attainable. 


1. Introduction 


The Fast Scattering Code (FSC) developed by Mark H. Dunn [ 1 1 requires the satisfaction of impenetrability condition on the 
scattering code. This condition requires the knowledge of the acoustic velocity produced by the incidence acoustic wave on 
the scattering surface. To find this velocity, one needs to find the gradient of the acoustic pressure on this surface. The 
particular incidence field that we are interested is generated by a helicopter rotor and is calculated from the rotor noise 
prediction code PSU-WOPWOP [2- 4], Any other code for noise prediction of rotating blades based on the acoustic analogy 
can also be used to generate the incident wave. 

Ideally, one would like to derive an analytic result for the gradient of the acoustic pressure which would require no more 
input data than those needed to calculate the thickness and loading noise. If we could get such a result, then the gradient of 
the acoustic pressure could be computed concurrently with the acoustic pressure using the geometric and kinematic data 
prepared for the acoustic calculations. This scheme will result in an efficient algorithm for calculating the gradient of the 
acoustic pressure on a computer. We will show in this paper that it is possible to obtain an analytic expression for the gradi- 
ent of the acoustic pressure which requires the same data as the thickness and loading noise for its evaluation. In the next 
section, we will present our derivation of this result. 


2. Derivation of The Gradient of the Acoustic pressure 


2.1. The Governing Wave Equations 

The governing wave equations that we will work with are those generating the thickness and loading noise in the Ffowcs 
Williams-Hawkings (FW-H) equation [5]. These are 
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Pt = 4~lPo v„ S(f)] (1) 

Ot 

□ 2 Pl=~ — iP n i 6(f)] ( 2 ) 

UX{ 

p = p T + p L (3) 

The blade surface is described by the equation f(x , t) = 0, where /(x, r) is defined in such a way that V/ = n, n being the 
unit outward normal to the blade with components n, . Here p is the acoustic pressure, and p T and p' L are the thickness and 
loading noise, respectively. The density of the undisturbed medium is p a and the local normal velocity of the blade surface is 
v„ . The blade surface pressure is denoted p and 6(f) is the Dirac delta function with the support on the blade surface / = 0 . 
The solution of the above wave equations was given by Farassat and Succi [6] and Brentner [7], The resulting solution, 
known as Formulation 1A of Farassat, is used in PSU-WOPWOP. The readers need some familiarity with the derivation of 
this Formulation to follow our mathematical work below. We recommend reference [7] for this purpose. 

2.2. The Gradient of the Acoustic Pressure 

We have 


V /V = V p T + V p L (4) 

We will now work with the formal solutions of eqs. (1) and (2) and manipulate them without integrating the delta functions 
in the integrands for several steps. We first work on finding the gradient of the thickness noise. 

Let (x, t) and (y, r) be the observer and the source space-time variables, respectively. The Green's function of the wave 
equation in the unbounded domain is: 

G (y, r; x,t)=^~ r<t 

^ A ” r (5) 

= 0 T > t 

g = T-t + - ( 6 ) 

c 

where r = \ x — y| and c is the speed of sound in the undisturbed medium. 

2.2.1. The Gradient of the Thickness Noise 

We introduce a new symbol for the thickness source strength as follows: Q = p a v„ . The formal solution of eq. (1) is: 

4 n p T (*’ t) = Jy <5(f) Kg) dr (7) 

where the volume integration is over the entire 3-dimensional space and the source time integration is over (-oo, /] . 

From eq. (7) we find that 

4 7tVp r (x,t) = V-^ J y 6(f) 6(g) dy dr = -^ f Q 6(f) V x (^dy dr ( 8 ) 

where the symbol V v stands for gradient operator with respect to the observer variable x. We now use the following result in 
the integrand of the last integral on the right of eq. (8): 

where r = (x - y) / r is the unit radiation vector from the source to the observer. Using this relation in eq. (8) and bringing the 
observer time derivative out of the integral, we obtain 
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<5/1 <5 CrQ 


£(/) 3(g) dy dr + 


6(f) 6(g) dydr^ 


( 10 ) 


The interpretation of the integrals in this equation using generalized function theory is fairly straight forward [5- 8]. We will 
just write down the final result here: 
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where M r is the local Mach number in radiation direction and 
Ei = 
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We now use the following result in the first integral on the right side of eq. (12): 
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Examining the integrand of the first integral on the right of the second equality sign, we note that 

Ej = 


1 d 

Q 'll 

1 -M r dr 

l r(l — M r ) )\ 


(15) 


is the integrand of the thickness noise [6, 7]. This we will write down shortly in detail. Furthermore, we can show that [7] 

dr c(M r r-M) c(l-M,.)r c(r - M) 

— = = + (16) 

dr r r r 

where M is the local Mach number of the blade surface based on the sound speed c. We remind the readers that the blade 
surface velocity used in M is described in the frame fixed to the undisturbed medium. 


Using eqs. (15) and (16) in eq. (14), we get 
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From references [6, 7], we have 
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( 18 ) 


where the dot on a symbol denotes the rate of change in the source time variable. We note that, by eq. (11), we have 
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4 n V p T (x, t) = - 
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dt 
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J ret 



(19) 


where E\ is given by eq. (17). We propose that the observer time derivative in eq. (19) be taken numerically in PSU-WOP - 
WOP. Then we see that for calculating the gradient of the thickness noise, we do not need any more data than is needed to 
calculate the thickness noise. This is the result we were looking for. 


2.2.2. The Gradient of the Loading Noise 

We will now find the gradient of the loading noise. We write eq. (2) in vector form as follows: 

a 2 p L = -V-[pnS(f)] (20) 

The Green's function solution of this equation is 

47tp L (x,t) = -V>J'~8(f)8(g)dydT = -J p8(f tfV x ^^-Jdydr (21) 

We now use eq. (9) in the integrand of the last integral and then take the observer time derivative outside the integral. The 
result is 
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From this result, we get 
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We remind the readers that at this stage of algebraic manipulations, the observer and the source space-time variables are 
independent because none of Dirac delta functions have been integrated. This approach makes it easy to interpret the differen- 
tiation operators. Had we been working with the integrated results, we would be dealing with heavy algebraic manipulations 
and the differentiation operators would require careful interpretation. 


To carry the analysis in eq. (23) further, we use the following two relations which can be derived easily 

If V 


tf Vv 


( r 8(g) 

= tf V x 

( r) 

- 

{ r 


\ rj 

[ ? 8(g) ) 

= H* V, 

r } 


r 2 J 


Tfrr, n - 2 cos 0 r cos Or d 

8 (g) = 5 <5(g) <%) 


c r 


c r dt 


tfr r , n - 3 cos 0 r cos Or d 

+ — S ( g ) = 1 8(g) — 5(g) 


c t 


c r 2 dt 


(24) 

(25) 


where cos 0 = if r , i.e., 0 is the angle between local outward normal to the blade surface and the radiation direction. Substitut- 
ing eqs. (24) and (25) in eq. (23) and gathering terms, we obtain 
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We now manipulate the first integral in the curly bracket on the right side of this equation to pull out the integrand of the 
loading noise. To help the readers to follow the algebraic steps below, we mention that the integrand of the loading noise is 
[6, 7]: 


E, = 


1 

1 5 1 

' p cos 9 t 

_L 

p cos 9 

c 

1 -M r dr l 

, r(l -M r ) ) 

-*ret 

_ r 2 (1 -M r ) 


p cos 9 


p ( cos 9 - M n ) 

1 

_l_ 

M» f p cos 9 

_l_ 

p cos 9 (M r - M 2 ) 

r ( 1 - M r ) 2 

T 

ret 

r 2 (1 - M,-) 2 

i 

c 

ret 

r (1 — M r ) 3 _ 

T 

ret 

r 2 (1 - M,-) 3 


(27) 


With this knowledge, we manipulate the following integral as follows 
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Here, we have used eq. (16) on the right of the first equality sign. We now use this result in eq. (26) which looks without 
simplification as: 
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We next simplify this result to get the result we are seeking: 
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We propose again that the observer time derivative in eq. (30) be taken numerically in PSU-WOPWOP. We see that we do 
not need any more data than those needed to calculate the loading noise to find the gradient of the loading noise. 


3. Concluding Remarks 


This analysis was performed to find the gradient of the acoustic pressure in PSU-WOPWOP which will be used in the Fast 
Scattering Code of Dunn [1]. Our main results, eqs. (19) and (30), use only the geometric and kinematic data available for 
thickness and loading noise calculation. This means that the gradient of the acoustic pressure can be calculated concurrently 
with the thickness and loading calculation in PSU-WOPWOP. The inclusion of our results in PSU-WOPWOP is quite easy 
and minimal programming will be involved. One does not expect noticeable increase in the execution time on a computer as 
compared to that of acoustic calculation only. 

We mention here that the choice of our main results was dictated by the use of Formulation 1A in PSU-WOPWOP. Another 
option available to us is to take the time derivatives in eqs. (19) and (30) analytically using the same method we have fol- 
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lowed above. Such a method will give a new analytic expression for the gradient of the acoustic pressure without the need for 
numerical differentiation with respect to the observer time. We have, however, rejected this procedure for the following 
reasons: 

1. The new analytic expression will have many more terms than those in eqs. (19) and (30). Checking the accuracy of the 
algebraic manipulations used in the derivation will be much harder. Programming and debugging of the new expression will 
be time consuming. 

2. There will be some new quantities needed in the evaluation of the new analytic expression, e.g., the second time derivative 
of the blade surface pressure, that will have to be prepared in the preprocessor. These quantities are not needed in Formula- 
tion 1 A and the preparation of these quantities for the input file will need new programming effort and debugging. 

3. Although numerical differentiation is avoided and this will reduce the execution time on a computer, more time will be 
spent on the evaluation of quantities that are not needed in Formulation 1A which will increase the execution time. There- 
fore, there may not be any net reduction of execution time by getting rid of the numerical differentiation. 

The extreme simplicity of implementing eqs. (19) and (30) in PSU-WOPWOP makes them ideal for the calculation of the 
gradient of the acoustic pressure needed in the Fast Scattering Code (FSC). 
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